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Abstract 

The frequency band used by cellular telephones includes the eigenfrequencies of a dielectric sphere 
with physical characteristics close to those of a human head. Proceeding from the spatial features 
of the natural modes of such a sphere we propose an independent and clear evident accuracy test 
for the complicated numerical calculations which are conducted when estimating the potential 
hazard due to the use of cellular telephones, in particular, for the check of a proper handling of 
the electromagnetic shape resonances of a human head. 
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Introduction. Estimation of the potential health hazard due to the use of cellular 
telephones is now a problem of primary importance in view of extremely rapid development 
and very wide spread of this communication aid. The safety guidelines in this field [l] 
are based on the findings from animal experiments that the biological hazards due to radio 
waves result mainly from the temperature rise in tissues 1 and a whole-body-averaged specific 
absorption rate (SAR) below 0.4 W/kg is not hazardous to human health. This corresponds 
to a limits on the output from the cellular phones (0.6 W at 900 MHz frequency band and 
0.27 W at 1.5 GHz frequency band). Obviously, the local absorption rate should be also 
considered especially in a human head [3]. Theoretical estimation of the temperature rise 
in tissues of a human body are accomplished in the following way. First the electric and 
magnetic fields inside the human body are calculated by solving the Maxwell equations with 
a given source (antenna of a portable telephone). The electric field gives rise to conduction 
currents with the energy dissipation rate aE 2 /2, where a is the conduction constant of 
respective tissues. In turn it leads to the temperature increase. The second step is the 
solution of the respective heat conduction equation (or more precisely, bio heat equation Q]) 
with local heat sources a E 2 /2 and with allowance for all the possible heat currents. Hence, 
for this method the distribution of electric field inside the human body (especially inside the 
head) is of primary importance. 

In this letter we are going to argue that the knowledge of the properties of electromagnetic 
modes for a dielectric sphere with physical characteristics close to those for a human head 
can be useful, for example, when developing an independent accuracy test of complicated 
numerical calculations mentioned above. The point is the eigenfrequencies of such a dielectric 
sphere lie in the GSM 400 MHz frequency band which has been used in a first generation 
of mobile phone systems and now is considered for further use. Obviously, the natural 
modes of a human head belong to this band too. The natural modes of a dielectric sphere 
can be divided into two types, surface and volume modes. For the volume modes the 
electromagnetic energy is distributed in the whole volume of the sphere while in the case 
of surface modes this energy is located close by the sphere surface. All this holds for the 
natural modes of a human head also, however we have no respective analytic formulas in 
this case. 



1 In principle, non- ionizing radiation can lead also to other effects in biological tissues Q. 



2 



In order to be fully confident, that the pertinent numerical schemes handle the resonances 
of a human head in a proper way, we propose an independent accuracy test of these calcu- 
lations. Without such a check it is not obvious because the routines of numerical solving 
the partial differential equations are local ones while the spatial behaviour of the relevant 
eigenf unctions characterizes the system whole. 

Shape resonances of a dielectric sphere. Let us consider a sphere of radius a, 
consisting of a material which is characterized by permittivity e\ and permeability The 
sphere is placed in an infinite medium with permittivity Ei and permeability yU 2 . It is assumed 
also that the electric currents are absent in both the media. The finite conductivity of the 
material inside a sphere will be taken into account below. 



It is known that in the source-free case (j = 0, p — 0) the general solution of Maxwe 
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equations are obtained from two scalar functions which may be chosen in different ways 
In the case of spherical symmetry these functions are the scalar Debye potentials ip (see, for 
example, the textbooks p, 13]): 

E™ = V x V x (rip™), H™ = -iuVx (rO (TM - modes), 

E™ = ico V x (r^™), H™ = V x V x (r^ E ) (TE - modes) . (1) 

The time dependence factor e~ lujt is dropped for simplicity. These potentials obey the 
Helmholtz equation inside and outside the sphere (r ^ a) and have the indicated angular 
dependence 

(V 2 + fc?M m = 0, fc f 2 = £i /i~, i = 1,2, il> lm (r) = f l (r)Y lm {n). (2) 

cr 

Equations (jSJ) should be supplemented by the boundary conditions at the origin, at the 
sphere surface and at infinity. In order for the fields to be finite at r = the Debye 
potentials should be regular there. Our goal is to find the eigenfrequencies and eigenfunctions 
in the problem at hand. Therefore at the spatial infinity, radiation conditions should be 
imposed js, i| 

lim r 1/2 fi{r) = const , lim r 1/2 ( - ik 2 fi(r)) = . (3) 

At the sphere surface the standard matching conditions for electric and magnetic fields 
should be satisfied jf]]. 
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In view of all this the Helmholtz equation becomes now the spectral problem for the 
Laplace operator multiplied by the discontinuous factor —l/(eiHi) 

1 oj 2 

&ipu>im(r) = -T^im(r), r^a, i= 1,2. (4) 



£{ fX{ C 

In this problem the spectral parameter is u 2 /c 2 . Due to the radiation conditions fl3l) this 

nn 

parameter is complex [9|, Ll0j - Thus we are dealing here with shape resonances of a dielectric 
sphere and the respective eigenf unctions are the quasi-normal modes 

In order to obey the boundary conditions at the origin and at spatial infinity formulated 
above, the solution to the spectral problem (jlj) should have the form 

f u i(r) = C l ji(k 1 r) , r<a, / wl (r) = C 2 hf\k 2 r) , r>a, (5) 

where ji(z) is the spherical Bessel function and h[ (z) is the spherical Hankel function of 



the first kind |l4j, the latter obeys the radiation conditions (02). 

At the sphere surface the tangential components of electric and magnetic fields ([1]) are 
continuous. As a result, the eigenfrequencies of electromagnetic field for this configuration 
are determined 6j by the frequency equation for the TE-modes 

y/ei^j'iiha) hi(k 2 a) - y/e 2 [J,iji(kia) ti^faa) = (6) 

and by the analogous equation for the TM-modes 

vfe 2 /ii j'i(ha) hi(k 2 a) - y/e^jiiha) h[(k 2 a) = , (7) 

where ki = y/Eifiiu/c, i = 1, 2 are the wave numbers inside and outside the sphere, 
respectively, and ji(z) and hi(z) are the Riccati-Bessel functions [yj] 

j l (z) = zji(z), h l {z) = zhf\z). (8) 

In equations and (J7J) the orbital momentum I assumes the values 1,2,..., and prime 
stands for the differentiation with respect of the arguments k±a and k 2 a of the Riccati- 
Bessel functions. 

The frequency equations for a dielectric sphere of permittivity e placed in vacuum follow 
from (jHJ) and (J7J) after putting there 

e x = e, e 2 = hi = fi 2 = 1 • (9) 



The roots of these equations have been studied in the Debye paper [15] by making use of an 
approximate method. As the starting solution the eigenfrequencies of a perfectly conducting 
sphere were used. In this case the frequencies are different for electromagnetic oscillations 
inside and outside sphere. Namely, inside sphere they are given by the roots of the equations 
(i>l) 

ji(^a) =0 (TE-modes), (10) 
^- (rji (^a)) = 0, r = a (TM-modes), (11) 

while outside sphere these frequencies are determined by equations 



h ( i 1) (^-aj=0 (TE-modes) (12) 
— (rh^ (^r)) = 0, r = a (TM-modes). (13) 

The frequency equations for perfectly conducting sphere ffTOl . ffTTj) and (JT2"]h fflBl can be 
formally derived by substituting ([9} into frequency equations © and © and taking there 
the limit e — > oo. 

Approximate calculation of the eigenfrequencies of a dielectric sphere without using com- 



puter 15| didnot allow one to reveal the characteristic features of the respective eigenfunc- 



tions (quasi- normal modes). The computer analysis of this spectral problem was accom- 
plished in the work where the experimental verification of the calculated frequencies 
was conducted also by making use of radio engineering measurements. 

These studies enable one to separate all the dielectric sphere modes into the interior and 
exterior modes and, at the same time, into the volume and surface modes. It is worth noting 
that all the eigenfrequencies are complex 

uj = uj' — i uj" . (14) 

Thus we are dealing with "leaky modes". It is not surprising because we are considering 
here an open system 1^] (a dielectric ball and outer unbounded space). 

The classification of the modes as the interior and exterior ones relies on the investigation 
of the behaviour of a given eigenfrequency in the limit e — > oo. The modes are called 
"interior" when the product ka = ^feujajc remains finite in the limit e —>■ oo, provided the 
imaginary part of the frequency (uj") tends to zero. The modes are referred to as "exterior" 



when the product kaj^fe = uoa/c remains finite with growing uj" . In the first case the 
frequency equations for a dielectric sphere (JBj) and (jTj) tend to ffTUj) and (TXTT) and in the 
second case they tend to ( 1121) and (fl3l) . The order of the root obtained will be denoted by 
the index r for interior modes and by r' for exterior modes. Thus TE/ r and TM; r denote 
the interior TE- and TM-modes, respectively, while TE; r / and TM; r / stand for the exterior 
TE- and TM-modes. 

For fixed I the number of the modes of exterior type is limited because the frequency 
equations for exterior oscillations of a perfect.y conducing sphere Q and « have finite 
number of solutions [9j. In view of this, the number of exterior TE- and TM-modes is given 
by the following rule. For even I there are 1/2 exterior TE- modes and 1/2 exterior TM- 
modes, for odd I the number of the modes TE/ r / is (/ + l)/2 and the number of the modes 
TM, r / equals (I- l)/2. 

An important parameter is the Q factor 

Q a - U ' -2-Ti stored energy 
ra 2uj" radiated energy per cycle 

For exterior modes the value of Qrad is always less than 1, hence these modes can never be 
observed as sharp resonances. At the same time for e greater than 5, the Q raa for interior 
modes is greater than 10 and it can reach very high values when e — >• oo. 

In the problem at hand the losses due to the radiation can be disregarded unlike the 
Ohmic losses. Indeed, the external source of electromagnetic energy (cellular telephone) 
compensates the radiation losses. While the Ohmic losses lead to the temperature rising in 
human tissues. 

For physical implications more important is the classification in terms of volume or surface 
modes according to whether r > I or I > r. For volume modes the electromagnetic energy 
is distributed in the whole volume of the sphere while in the case of surface modes the 
energy is located close by the sphere surface. The exterior modes are the first roots of the 
characteristic equations and it can be shown that they are always surface modes. 

Figure 1 shows a typical spatial behaviour of the surface and volume modes of a dielectric 
sphere. 

Thus a substantial part of the sphere modes (about one half) belong to the interior surface 
modes. It is important that respective frequencies are the first roots of the characteristic 
equations. 



6 




0.0 0.2 0.4 0.6 




0.0 0.2 



FIG. 1: Electric energy density r 2 E 2 for the surface (A) and volume (B) TE-modes of a dielectric 
sphere with e = 40 placed in vacuum. 

In order to escape the confusion, it is worth noting here that the surface modes in the 
problem in question obey the same boundary conditions at the sphere surface and when 
r — > oo as the volume modes do. Hence, these surface modes cannot be classified as the 
evanescent surface waves propagating along the interface between two media (propagating 
waves along dielectric wave guid es 



bulk and adjacent vacuum 
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7| , surface plasmon waves on the interface between metal 
19| and so on). When describing the evanescent waves one 
imposes the requirement of their exponential decaying away from interface between two 
media. In this respect the evanescent surface wave differ from the modes in the bulk. 

Features of dielectric sphere spectrum and their applications. The parts of 
human body (for example, head) have the eigenfrequencies of electromagnetic oscillations 
like any compact body. In particular, one can anticipate that the eigenfrequencies of human 
head are close to those of a dielectric sphere with radius a ~ 8 cm and permittivity e ~ 40 
(for human brain e = 44.1 for 900 MHz and e = 42.8 for 1.5 GHz [3]). By making use of the 
results of calculations conducted in the work 16] one can easily obtain the eigenfrequencies 
of a dielectric sphere with the parameters mentioned above. For TEji modes with I = 1, 2, 3 
we have, respectively, the following frequencies: 280 MHz, 420 MHz, and 545 MHz. For TMn 
modes with / = 1, 2, 3 the resonance frequencies are 425 MHz, 540 MHz, and 665 MHz. The 
imaginary parts of these eigenfrequencies are very small so the Q factor in f|T5|) responsible 
for radiation is greater than 100. 

These eigenfrequencies belong to a new GSM 400 MHz frequency band which is now 
being standardized by the European Telecommunications Standards Institute. This band 
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was primarily used in Nordic countries, Eastern Europe, and Russia in a first generation of 
mobile phone system prior to the introduction of GSM. 

Due to the Ohmic losses the resonances of a dielectric sphere in question are in fact 
broad and overlapping. Indeed, the electric conductance a of the human brain is rather 
substantial. According to the data presented in Q| a ~ 1.0 S/m. The eigenfrequencies of 
a dielectric dissipative sphere with allowance for a finite conductance a can be found in 



the following way. As known 20] the effects of a on electromagnetic processes in a media 
possessing a common real dielectric constant e are described by a complex dielectric constant 
£diss depending on frequency 

Ana 

£diss = £ + i • (16) 

The eigenfrequencies u>, calculated for a real e, are related to eigenfrequencies a;di ss for e^ss 



by the formula 
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\ 1/2 

^diss = ( I uj^uj~2iri-. (17) 

Sdiss/ e 



The corresponding factor Qdiss is 

2 ^diss 47r 

Substituting in this equation the values u>/2n = 0.5 • 10 9 Hz, e = 40, a = lS/m = 
9 ■ 10 9 s _1 one finds 

20 

Qdte- — ~1. (19) 

In view of such substantial Ohmic loses the resonance enhancement of the oscillation 
amplitude inside a human head will not occur. However when the frequency of a mobile 
telephone coincides with the eigenfrequency of the head the distribution of electric and 
magnetic fields inside the head will be described by the corresponding normal mode which 
may be a surface mode or a volume one 
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Proceeding from this we propose the following test of numerical calculations used when 
estimating the potential hazard of cellular telephones. The test consists in simulation of the 
temperature distribution corresponding to the surface and volume modes in the framework 
of pertinent calculation schemes. For simplicity, the test calculations could be accomplished 
for a dielectric sphere (instead of a human head) with lower conductivity in comparison 
with that for a human brain (in order to enhance the effect). The distributions of electric 
and magnetic fields and the temperature distribution inside the sphere should be calculated 
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for two eigenfrequencies of the sphere, namely, one frequency corresponds to surface mode 
and another one belongs to volume mode. The distributions obtained should conform, at 
least qualitatively, to the spatial behaviour of respective electromagnetic normal modes (see 
Fig. 1). 

Conclusion. Detailed analysis of electromagnetic spectra of a dielectric sphere enables 
us to propose an independent accuracy test of complicated numerical calculations conducted 
when estimating the potential health hazard due to use of cellular telephones. This test will 
permit one to make certain of a proper handling of the electromagnetic shape resonances of 
a human head in these studies. 

This paper was completed during the visit of on of the authors (VVN) to Salerno Univer- 
sity. It is his pleasant duty to thank G. Scarpetta and G. Lambiase for the kind hospitality 
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